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SUMMARY 

1  1 

The  simulation  studies  described  shott  that  sign-logarithm  arithmetic  can  be 
implemented  in  a  practical  digital  FfT  analyser.  Sign-logarithm  arithmetic 
allows  a  smaller  wordlength  than  conventional  fixed  point  arithmetic  whilst 
maintaining  performance. 


Discussion  of  the  hardware  implementation  of  such  a  sign-logarithm  FFT  shows 
that  power  consumption  can  be  less  than  conventional  methods  using  bipolar 
multipliers.  The  use  of  a  smaller  wordlength  allows  a  significant  simplifi¬ 
cation  of  the  system  into  which  the  FFT  is  placed  and  a  higher  data  throughput 
rate. 
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ABSTRACT 

The  simulation  studies  described  show  that  sign-logarithm  arithmetic  can 
be  implemented  in  a  practical  digital  FFT  analyser.  Sign  logarithm  arithmetic 
allows  a  smaller  wordlength  than  conventional  fixed  point  arithmetic  whilst 
maintaining  performance. 

Discussion  of  the  hardware  implementation  of  such  a  sign-logarithm  FFT 
shows  that  power  consumption  can  be  less  than  conventional  methods  using  bi¬ 
polar  multipliers.  The  use  of  a  smaller  wordlength  allows  a  significant 
simplification  of  the  system  into  which  the  FFT  analyser  is  placed  and  a  higher 
data  throughput  rate. 
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1  INTRODUCTION 

Traditional  hardware  implementations  of  FFT  algorithms,  using  fixed  point 
arithmetic,  suffer  from  dynamic  range  limitations  and  slow  multiply  speeds. 

The  problem  of  multiply  speed  becomes  more  acute  with  the  longer  wordlengths 
required  for  a  good  dynamic  range. 

Sign-logarithm  arithmetic  attempts  to  overcome  these  problems  by: 

1  Increasing  the  dynamic  range  for  a  given  number  of  bits  by  non-linear 
quantisation. 

2  Reducing  arithmetical  complexity  and  time,  since  multiplication  can  be 
replaced  by  addition  of  base  2  logarithms. 
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The  scope  of  this  memorandum  is  to  show  how  sign-logarithm  arithmetic 
functions  and  how  it  can  be  implemented  in  a  ITT  analyser.  The  sign-logarithm 
FFT  was  simulated  and  compared  with  a  fixed  point  FFT.  A  feasibility  study  was 
then  conducted  to  determine  the  usefulness  of  this  technique  in  a  practical 
system. 

2  GENERALISED  LOGARITHM  ARITHMETIC 
2.1  MULTIPLICATION 

The  expression  used  for  multiplication  by  addition  of  logarithms  is 
explained  in  most  standard  algebra  texts. 


let 


z  -  x  *  y 

a  =  l°gp  W 

b  «  logp|y| 

Sa  «  sign(x) 

where  S^  «  1,  k 

Sb  =  sign(y) 

-  -1,  1 

-  0,  k 

m  -  logp|z| 

"  logp(|x  *  y|) 

Sm  *  sign(z) 

«  sign(x  *  y) 

then 


2.2 


m  -  logpjx|  +  logp |y | 
Sm  «  sign(x)  *  sign(y) 

ADDITION  AND  SUBTRACTION 


Addition  and  subtraction  can  be  achieved  when  the  logarithm  of  the  two 
numbers  is  known,  and  an  expression  is  easily  obtained:  [1] 


let  s 
a 
b 
Sa 
Sb 

then  c 
Sc 


x  +  y 
logp  U I 
logplyl 
sign(x) 
sign(y) 

logpM  * 
sign(s)  « 

(i)  Assume  that  Sa 

then  Sa 


logp (|x  +  y!) 

8ign(x  +  y) 

■  Sb. 

■  Sb  *  Sc 

■  log  (p*  +  Pb> 

r  b— * 

■  a  +  logpCl  p  ) 

■  b  ♦  logp(l  ♦  pa  k) 


2 


letting  d  -  | a  -  b] 

then  IF  a  >  b,  c  ■ 
IF  a  <  b,  c  ■ 
IF  a  -  b,  c  - 

m 

The  correction  factor, 
set  of  values;  since: 


a  +  logp(l  +  p  d) 
b  +  logp(l  +  p~d) 
a  +  logp(2) 

b  +  logp (2) 

dog  d  +  p  d))  can 


only  have  a  limited 


0  <  p"d  <  1 


0  <  log  Cl  +  p~d)  <  log  2 
P  P 

(ii)  Assume  that  Sa  «  -Sb  and  a  >  b 


Using  the  same  notation  as  before: 
since  a  >  b  then  |x|  >  |y| 

c  -  logp(|x  +  y|) 

“  logpdxl  -  | y t ) 

-  log  (pa  -  pb) 

P  -d 

-  a  +  logpd  -  P  ) 

Similarly,  if  b  >  a  then: 
c  »  b  +  logp(l  -  p  d) 

If,  however,  b  »  a,  then  the  correct  answer  is  S  «  0,  and  a 
special  case  has  to  be  made  for  the  value  of: 

c  -  logp(O) 

Addition  and  subtraction  can  be  accomodated  by  adding  a  correc¬ 
tion  factor  to  the  largest  operand.  The  correction  factor  (6)  is 
calculated  by  the  functions: 

{  -  B(a,b)  ■>  logpd  +  P_d)  if  Sa  -  Sb 

-  y(a,b)  •  logpd  -  p  d)  if  Sa  -  -Sb 

where  d  ■  ja  -  b| 

3  FIXED  POINT,  BINARY  SIGN-LOGARITHM  ARITHMETIC 

The  methods  developed  in  the  previous  section  need  modification  to  accommo¬ 
date  fixed  wordlengths.  In  general,  a  binary  word  can  be  split  into: 


1  sign  bit 
p  integer  bits 
q  nominal  fraction  bits 


for  example  011.010  could  represent  3.25 
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In  the  sign-logarithm  number  system,  a  real  number  a  is  represented  by  its 
sign  (Sa)  and  the  logarithm  of  its  magnitude  La. 

IF  la1  >  I  THEN  h&  -  log2(x|a|) 

IF  a  <  -  THEN  L  -  0 
x  a 

IF  a  <  0  THEN  S  -  1 
a 

IF  a  >  0  THEN  S  -  0 
a 

If  q  bits  are  used  to  represent  the  fractional  part  of  the  logarithm  of  a, 
La  may  be  expressed  thus: 

IF  la  I  >  ~  THEN  Lg  -  ENTIERly  +  2qlog2  (x  |a  I )]  /2q 

The  scale  factor  x  ensures  positive  logarithms  for  small  numbers.  The 
constant  of  1/2  unbiases  the  quantisation  error  by  rounding,  rather  than 
truncating. 

Arithmetical  calculations  can  be  handled  as  before: 


3. 1  MULTIPLICATION 

P  -  a  *  b  Lg  %  log2(x  ja() 

1^  ■»  log2(x  |b|) 

Lp  %  1o82(t'p^ 

Lp  ’  lo82^a  *  bl) 

»  log2(x|a|  *  x  |b  |  *  i) 

-  log2(x|a|)  +  log2 (x |b  |)  -  log2(x) 

«  La  +  Lb  -  Lc 

where  Lc  is  a  correction  factor  Lc  -  ENTIERfy  +  2qlog2 (x)J  /2q 
The  sign  of  the  product  is  determined  by: 

Sp  *  Sa  ®  Sb 

The  expression  Lp  »  La  +  Lb  -  Lc,  allows  both  underflow  and  overflow 
to  occur. 


Underflow  occurs  when  Lp  is  negative,  and  indicates  that  the  correct 
product  is  smaller  than  the  smallest  sign-logarithm  number,  ie  p  <  1/x. 

It  is  useful  to  detect  this  underflow,  and  assign  zero  to  the  result. 

Overflow  occurs  when  the  result  exceeds  the  longest  number  that  can 
be  represented  by  an  n  bit  binary  logarithm. 

ie  Lp  >  2n_q  -  1 

Numbers  can  be  scaled  by  right  shifts,  and  doubling  x;  or  subtracting 
a  constant  from  the  logarithm. 
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3.2  ADDITION  AND  SUBTRACTION 
s  *  a  +  b 

(i)  Assuming  Sa  «  Sb  and  Lfl  > 

from  the  generalised  analysis  c  «  log  (la  +  b|) 

p  -d 

-  a  +  logp(l  +  p  ) 

80  Ls  “  La  +  6avV 
where  6(La,L^)  is  defined  by 

6(^,1^)  -  ENTIERf-i  +  2qlog2(l  +  2'd» /2q 

d  ‘  'La  '  S' 

The  function  8(La,L^)  has  a  limited  set  of  values  in  the  region 

(0,1). 

(ii)  Assuming  Sa  ■  -Sb  and  a  >  b 
La>Lb 

From  the  generalised  analysis  c  «  log  (|a  +  b|) 

»  a  +  logp(l  “  P  d> 

80  Ls  "  La  *  y(W 
where  >(^,1^)  is  defined  by 

Yaa,l1j)  -  ENTIERlj  +  2qlog2(l  -  2_d)]  /2q 
d  -  |La  -  Lj 

Since  d  >  0,  a  special  case  exists  when  d  «  0,  ■  L^. 

IF  d  -  0  THEN  Ls  -  0 

The  sign  of  the  solution  is  given  by  the  largest  value,  in  this 
case  Ss  ”  Sa. 

4  LOGARITHMIC  QUANTISATION 

Non-linear  quantisation  is  frequently  used  in  speech  communication  systems 
to  improve  dynamic  range  and  to  maintain  a  constant  signal  to  quantisation 
noise  power  ratio. 

Any  system  of  sign-logarithm  arithmetic  requires  quantisation  of  the 
signal  in  logarithmic  intervals.  Thus  the  advantages  of  non-linear  quanti¬ 
sation  also  apply  to  any  system  using  sign-logarithm  arithmetic. 

A  first  order  approximation  for  the  quantisation  noise  for  a  sinusoid  is 
as  follows. 
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Substituting  for  dx  where  x  »  t  and  large  m  gives 
1 


Nq 


12  (to) 
1 

3K2m2  f. 


p(x)dx 


x  p(x)dx 


3K  m 


S  o23K2m2 
Thus  ^  »  - j- 


since  we  have  assumed  that  th 
wholly  contained  in  the  inten 


3K2m2 


smal  is 


for  x  >  t 

and  is  independent  of  the  signal. 


A  more  detailed  analysis,  which  accounts  for  all  signal  levels,  is 
shown  in  Figure  2  (with  l  =  6  x  lO^  for  the  logarithm  quantisation).  It 
can  be  seen  that  the  signal  to  quantisation  noise  power  for  logarithmic 
quantisation  remains  constant  over  a  much  wider  range  than  linear  quanti¬ 
sation. 


It  may  therefore  be  assumed  that  for  a  system  which  encounters  a  wide 
variety  of  signal  levels,  logarithm  quantisation  should  be  used. 

5  SIMULATION 

5.1  METHOD  OF  ANALYSIS 

A  sign-logarithm  Radix-2  FFT  processor  was  simulated  on  a  computer  to 
verify  this  technique.  The  simulation  was  based  upon  the  rules  described 
in  previous  sections.  The  sign-logarithm  processing  was  compared  with  a 
simulated  fixed  point  Radix-2  FFT,  for  different  wordlengths.  Both  FFTs 
automatically  rescaled  their  operands  when  it  was  likely  that  an  operation 
would  result  in  an  overflow.  Hamming  windowing  was  included  in  the  pro¬ 
cessing  of  both  FFTs. 

The  output  from  the  sign-logarithm  FFT  was  expanded  so  that  a  reason¬ 
able  comparison  could  be  made  (compression  removed).  The  ADC  clipping 
levels  were  set  at  ±1  V.  The  peak  signal  amplitude  was  0.25  V,  with  a 
signal  to  noise  ratio  of  40  dB. 

Results  were  displayed  graphically. 

5.2  SIMULATION  RESULTS 

The  results  are  described  briefly  in  this  section,  but  are  shown  in 
store  detail  in  Appendix  A. 
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(i)  For  a  fixed  wordlength,  the  number  of  bits  used  for  the  frac¬ 
tional  part  of  the  logarithm  affects  the  sidelobe  levels  of  the 
sign- logarithm  FFT.  An  optimum  number  of  fractional  bits  exists 
for  each  wordlength. 

WORDLENGTH  OPTIMUM  NUMBER  OF  FRACTIONAL  BITS' 

8  3 

6  2 

5  1 

(ii)  For  the  wordlengths  used,  and  with  an  optimum  number  of  frac¬ 
tional  bits,  the  sign-logarithm  FFT  gave  superior  results  to  the 
fixed  point  FFT  with  the  same  wordlength.  This  may  be  attributed 
to  the  increased  dynamic  range. 

5.3  DISCUSSION  OF  RESULTS 

From  the  simulation,  two  implications  exist: 

(i)  For  the  same  wordlengths,  superior  results  may  be  obtained  using 
a  sign-logarithm  FFT. 

(ii)  The  same  performance  as  a  fixed  point  FFT  may  be  achieved  with  a 
smaller  wordlength,  by  using  sign-logarithm  arithmetic. 

The  ability  to  add  logarithms  in  place  of  multiplication,  coupled 
with  the  improved  performance  of  the  sign-logarithm  FFT,  indicates  that 
this  method  could  offer  significant  advantages  in  processing  time  and 
power  dissipation  over  conventional  methods 

6  HARDWARE  IMPLEMENTATION 
6.1  ARITHMETIC  UNITS 

The  radix-2  butterfly  structure  requires  four  real  multipliers  and 
six  add/subtract  units  (Figure  3).  A  radix-2  butterfly  for  sign-logarithm 
arithmetic  may  be  constructed  from  hardware  which  follows  the  arithmetical 
rules  established  in  section  3,  assuming  that  the  data  has  been  quantised 
in  a  logarithmic  fashion  as  described  in  section  4. 

A  multiplier  requires  the  addition  of  the  logarithms  according  to  the 

rules 


Lp  “  La  +  Lfc  -  Lc 

Sp  -  Sa  ©  Sb  (3.1) 

The  required  hardware  i6  shown  in  Figure  4,  which  also  deals  with 
underflow  for  small  numbers.  Since  one  operand  is  a  weighting  function 
stored  in  ROM,  the  weighting  function  may  be  adjusted  prior  to  addition  so 
that  it  includes  the  correction  factor  Lc. 

An  add/subtract  unit  requires  a  more  complex  arrangement  shown  in 
Figure  5.  The  correction  factors  y(La,Lb)  and  B(La,Lb)  are  modified  to: 
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Y (d) ,  B(d) 

enabling  the  coefficients  to  be  stored  in  a  small  ROM  and  addressed  by  the 
difference  between  La  and  Lj,. 

6.2  RADIX-2  BUTTERFLY 

A  radix-2  butterfly  may  be  constructed  from  the  units  described  in  6.1 
as  shown  in  Figure  6.  Some  redundancy  has  enabled  the  component  count  to 
be  reduced.  Scaling  of  the  results  may  be  achieved  by  adjusting  the  look 
up  tables  in  the  final  add/subtract  units  as  shown. 

6.3  REQUIREMENTS  AND  CONSEQUENCES  OF  IMPLEMENTATION 

6.3.1  COMPARISON  OF  SPEED  AND  POWER 

Comparisons  were  made  on  the  basis  that  sixteen  bit  fixed  point 
FFTs  are  equivalent  to  eight  bit  sign-logarithm  FFTs.  The  figures 
shown  are  for  the  parallel  implementation  of  one  butterfly  shown  in 
Figure  3. 

a  Speed 

The  longest  propagation  delay  in  the  radix-2  butterfly  is 
the  tije  required  to  complete  one  multiply  and  two  adds.  The 
sign-logarithm  implementation  required  five  adds  and  two  ROM 
delays  (Figure  6) . 


bipolar 

16  bit 

100 

+  26 

♦ 

26 

* 

152 

ns 

cmos 

16  bit 

(TMC  216H) 

170 

+  26 

+ 

26 

m 

222 

ns 

cmos 

16  bit 

(ADSP  1016) 

130 

+  26 

♦ 

26 

m 

182 

ns 

bipolar 

16  bit 

(MPY  01 6K) 

45 

+  26 

+ 

26 

w 

97 

ns 

bipolar 

8  bit 

sign-logarithm 

5  x 

20  + 

2 

x  30 

m 

160 

ns 

b  Power 


Calculations  are  based  on  the  typical  power  consumptions 
for  all  the  components  required  to  implement  one  butterfly. 

This  requires: 

conventional  implementation  -  4  multipliers,  6  x  16  bit 

adders. 

sign-logarithm  implementation  -  14  x  8  bit  adders,  4  x  128 

x  8  ROMs. 

4x8  bit  MUX,  4x8  bit 
comparators. 


bipolar 

16  bit 

30  W 

cmos 

16  bit 

(TMC  216H) 

14  W 

cmos 

16  bit 

(ADSP  1016) 

12.6  V 

bipolar 

16  bit 

(MPY  016K) 

30  W 

bipolar 

8  bit 

sign-logarithm 

16  W 
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X  -  R  +  B*EXP(j2*PI*n/N> 
Y  -  R  -  B*EXP(j2#PI*n/'N) 


Figure  3.  Radix-2  Butterfly,  real  and  Imaginary  parts. 


.Figure  4.  Hardware  Implementation  of 
Sign-logarithm  multiplier. 
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6.3.2  DISCUSSION 
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The  gain  in  processing  time  made  by  adding  logarithms  for  multi¬ 
plication  is  offset  by  the  more  complicated  procedure  required  to  add 
two  numbers  together.  Whilst  the  sign-logarithm  butterfly  consumes 
less  power  than  its  equivalent  bipolar  fixed  point  butterfly,  this 
advantage  has  been  'leapfrogged'  by  fast  cmos  multipliers. 

The  sign-logarithm  FFT  is  likely  to  be  cheaper  to  implement 
whilst  the  cost  of  multiply  chips  remains  high. 

The  logarithmic  quantisation  may  be  achieved  by  preceding  an  ADC 
with  a  non-linear  amplifier  with  the  required  transfer  characteristic. 

Consideration  of  a  complete  system  into  which  a  sign-logarithm 
FFT  is  placed  shows  several  advantages. 

a  ADCs  with  fewer  bits  can  be  used.  These  are  faster  and 
consume  less  power. 

b  Smaller  wordlengths  generally  allow  faster  data  throughput 
rates  in  all  circuits. 

c  Smaller  wordlength  allows  simpler  board  layouts,  but 
additions  increase  the  complexity. 

d  The  greater  dynamic  range  available  with  non-linear  quanti¬ 
sation  allows  earlier  analogue  to  digital  conversion  in  a 
system.  The  need  for  filtering  to  reduce  dynamic  range 
before  conversion  is  eliminated. 

7  CONCLUSIONS 


The  sign-logarithm  FFT  is  an  attractive  alternative  to  conventional  fixed 
point  implementations.  It  offers: 


a  greater  dynamic  range  for  an  identical  number  of  bits 

b  almost  constant  signal  to  quantisation  noise  power  ratio 

c  a  reduction  in  power  consumption  when  coopered  to  bipolar  multiplier 
configurations 


d  1 ler  wordlengths  for  similar  performance 

f  cication  of  preceding  and  subsequent  signal  processing  because 

laller  wordlengths 


Thi 

radix-2  dull. 
add  two  numbers 


eved  with  a  slight  increase  in  the  time  required  to  perform  a 
,  a  consequence  of  the  more  complicated  procedure  required  to 
together. 
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APPENDIX  A 
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SIMULATION  RESULTS 

The  results  of  computer  simulation  of  fixed  point  FFTs  and  sign-logarithm 
FFTs  with  different  vordlengths  are  shown  in  Figures  A1  to  A28. 

The  figures  are  shown  in  two  formats. 

a  An  annotated  contour  plot  of  the  modulus  output  of  15  x  16  point  com¬ 
plex  FFTs.  Points  exist  on  the  plus  marks,  contours  are  interpolated 
between  them  at  -3  dB  intervals.  Plots  are  shown  for  many  different 
word  configurations. 

b  A  graph  of  the  modulus  output  of  one  or  more  16  point  FFTs  for  differ¬ 
ing  word  configurations. 

KEY  TO  FIGURES 

Al  Contour  plot  with  16  bit  fixed  point  arithmetic. 

A2  Contour  plot  with  8  bit  fixed  point  arithmetic. 

Not  significantly  different  to  Al  down  to  -30  dB. 

A3)  Contour  plot  with  8  bit  sign-logarithm  arithmetic. 

A4;  Optimum  number  of  bits  for  the  fractional  part  is  3 
A5J  and  shows  no  significant  difference  to  A2  down  to  -30  dB. 

A6  Contour  plot  with  6  bit  fixed  point  arithmetic. 

Sidelobe  levels  higher  than  -24  dB  due  to  the 
inability  to  resolve  the  full  dynamic  range  required. 

A7  )  Contour  plot  with  6  bit  sign- logarithm  arithmetic. 

A8>  Optimum  number  of  fractional  bits  is  2. 

A9 J  A  significant  improvement  in  sidelobe  levels  over  A6. 

A10  5  bit  fixed  point  arithmetic. 

Sidelobe  levels  higher  than  -18  dB. 

All)  5  bit  sign- logarithm  arithmetic. 

A12>  Showing  optimum  number  of  fractional  bits  is  1. 

A13;  A  significant  improvement  in  sidelobe  levels  over  A10. 

A14)  5  bit  arithmetic,  showing  an  18  dB  improvement  in 
A159  dynamic  range  available  with  sign-logarithm  arithmetic, 
with  no  loss  of  precision. 

A16)  Showing  a  significant  improvement  in  dynamic  range  using 
A17/  8  bit  sign- logarithm  arithmetic  over  8  bit  fixed  point. 

A18;  Performance  of  the  8  bit  sign-logarithm  FFT  is  comparable 
to  the  16  bit  fixed  point  FFT. 
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16  bit  fixed  point  vs  8  bit  sign-logarithm 


A19 

A20 

A21 

A22 

A23 


Contour  plot  and  graphs  of  the  modulus  output  of 
16  bit  fixed  point  ITT  with  more  than  one  frequency 
component . 


A24 

A25 

A26 

A27 

A28 


Contour  plot  and  graphs  of  the  modulus  output  of 
8  bit  sign-logarithm  FT!  with  more  than  one 
frequency  component. 


The  two  systems  do  not  differ  signif icantly  in  performance. 
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